Insights into the electronic structure of Co 2 FeSi from x-ray magnetic linear dichroism 
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Experimental evidence both for and against a half-metallic ground-state of the Heusler com- 
pound Co2FeSi has been published. Density functional theory based calculations suggest a non 
half-metallic ground state. It has been argued, that on-site Coulomb interaction of the d electrons 
has to be taken into account via the LDA+C/ method, which predicts a half-metallic ground-state 
for U ~ 2.5 . . . 4.5 eV. X-ray magnetic linear dichroism (XMLD) can be used as a tool to assess the 
appropriateness of the LDA+t/ approach: the calculated spectra within the LDA+C/ or GGA+C/ 
schemes are different from those within the LDA or GGA. Due to its ability to separate differ- 
ent orbital symmetries, XMLD allows us to distinguish between different models of the electronic 
structure of Co2FeSi. In this article we discuss experimental XMLD spectra and compare them with 
detailed first principles calculations. Our findings give evidence for the inadequacy of the LDA+(7 or 
GGA+C/ band structures, whereas constrained calculations with the GGA and a fixed spin moment 
of 6 flu give better overall agreement between experiment and theory. 



I. INTRODUCTION 

Half-metallic Heusler compounds of the type C02YZ 
(with Ya3rf transition metal and Z an sp main group el- 
ement) are promising materials for spin electronics. They 
were shown to follow the Slater-Pauling rule, which gives 
their magnetization per formula unit as m — N v — 24, 
where N v is the number of valence electrons^ Co2FeSi is 
the Heusler compound with the highest observed magne- 
tization of nearly 6/xb per formula uniti^ 

The electronic structures of half-metallic Heusler com- 
pounds are routinely computed within the framework of 
density functional theory (DFT). The local (spin) density 
and generalized gradient approximations (LDA or GGA) 
are often able to reproduce the observed ground-state 
magnetic moments in agreement with experimental data. 
However, Co2FeSi is a case, for which these approxima- 
tions fail: the calculated magnetic moment is too low and 
the density of states is not half- metallic. Because of this, 
approaches to go beyond the (semi-)local approximations 
were tested. Applying LDA+C/ in the fully localized limit 
with U around 4 eV opens a minority gap around E-p and 
consequently gives the full magnetic moment 3 Calcula- 
tions with the non-empirical hybrid PBEO functional re- 
sulted in a half-metallic ground state with a 1 eV wide 
minority gap. 3 The LDA+C/ (or GGA+C/) approach has 
been applied to study the electronic structure of Co2FeSi 
and related compounds 4^— 

Although Co2FeSi satisfies the Slater-Pauling rule and 
has a near- integer magnetic moment per formula unit, 
which is a necessary condition for half-metallic ferro- 
magnetism, there is some experimental evidence to the 
contrary. Tunnel magnetoresistance experiments show 
rather low TMR with Co2FeSi electrodes, and the con- 
ductance curves lack characteristic features of the mi- 
nority gap^iS Further, the Gilbert damping is quite 



highlit On the other hand, a temperature independent 
resistivity is observed for bulk single- and polycrystallinc 
Co2FeSi below 50 K, which is seen as a fingerprint of half- 
metallicitypii 

Recently, Kallmayer et al. predicted from first- 
principles calculations the x-ray magnetic linear dichro- 
ism (XMLD) of Co 2 FeSi within the LDA and the 
LDA+C/ schemes*^ No significant differences of the x-ray 
absorption or circular dichroism spectra were predicted. 
However, clearly different XMLD spectra were computed, 
and XMLD was proposed as a tool to determine which 
electronic structure is actually present in Co2FeSi. In 
this article we use the proposed method on epitaxial thin 
films of Co 2 FeSi. 

The article is organized as follows: Sections |H] and 
IIIII describe the experimental conditions and details of 
our computational approach. The results of the x-ray 
magnetic circular (XMCD) and linear dichroism mea- 
surements and calculations are presented in Section IIV1 
followed by a discussion of the results in Section [V] The 
main part of this article is focused on the understanding 
of the XMLD of Co2FeSi. To make sure that the dichroic 
signals are intrinsical to Co2FeSi, we investigated samples 
prepared with different deposition conditions. 



II. EXPERIMENT 

Epitaxial thin films of Co2FeSi have been grown 
on MgO(OOl) substrates by dc and rf magnetron co- 
sputtering with and without seed layers and various heat 
treatments. 

The samples were deposited with a co-sputtering tool 
with 3" sputter sources. It allows heating of the 
substrates to 900° C. The base pressure was typically 
2 ■ 10~ 9 mbar, and the Ar working pressure was 2 • 
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TABLE I. Preparation conditions of the samples. 



Name Stack (thicknesses in nm) 



CFS-Cr MgO(OOl) / Cr 5 / 700°C in situ / CFS 20 @ RT/ Mg 0.5 / MgO 1.0 / 450°C ex situ 

CFS-HTD MgO(OOl) / CFS 20 @ 700°C / Mg 0.5 / MgO 1.0 



lCP 3 mbar. The target-to-substrate distance was 21cm. 
Elemental Co, Fe, and Si targets of 99.9% purity were 
used. The metals were dc sputtered and Si was rf sput- 
tered. The MgO capping layers were deposited by elec- 
tron beam evaporation. 

We distinguish two types of samples, which are ex- 
plained in Table |U CFS-Cr (with Cr seed layer and post- 
annealing) and CFS-HTD (high-temperature deposition 
without seed layer). 

X-ray absorption (XAS), x-ray magnetic circular 
(XMCD) and linear dichroism (XMLD) spectra were 
taken at beamline 4.0.2 of the Advanced Light Source, 
Berkeley. Magnetic fields of 0.55 T were applied to mag- 
netically saturate the samples. XMCD spectra were 
taken with an angle of 30° between the beam and the 
sample surface, the magnetic field was parallel to the 
beam. The circular polarization degree was 90%. XMLD 
spectra were taken in normal incidence, with the mag- 
netic field canted out of the sample plane by 10°. The 
degree of linear polarization was 100%. The magnetiza- 
tion was switched at every photon energy to obtain the 
dichroic spectra. All measurements were taken at room 
temperature, which is justified in view of the high Curie 
temperature of Co2FeSi. Surface sensitive total electron 
yield spectra (TEY) were taken concurrently with trans- 
mission spectra. The transmission signal was collected 
by a photodiode behind the sample, which measured the 
x-ray luminescence in the substrate The beamline res- 
olution was about 0.1 eV. 

Structural characterization was performed by x-ray 
diffraction in a Philips X'Pert Pro MPD diffractometer 
with a Cu anode. Bragg-Brentano optics were used for 
the specular measurements, and an open Eulerian cradle 
in combination with collimator point-focus optics was ap- 
plied for off-specular measurements. 

Film compositions were determined by x-ray fluo- 
rescence spectroscopy, which has been cross-checked 
with inductively coupled plasma optical emission spec- 
troscopy. The film stoichiometries were thus adjusted to 
Co5o±2Fe 2 5±iSi 2 5±i- 

Magnetization measurements have been performed 
with a vibrating sample magnetometer at 50 K and 300 K 
in magnetic fields up to 3 T. 



III. COMPUTATIONAL APPROACH 

The computational work was carried out within the 
Elk code^ which is based on the full potential linearized 
augmented plane waves (FLAPW) method. Various ap- 



proximations to the exchange-correlation functional were 
tested. Within the generalized gradient approximation 
(GGA), we chose the Perdew-Burke-Ernzerhof (PBE) 
functional^ We performed LDA+J7 calculations, based 
on the Perdew-Wang local (spin) density approximation 
(LDA). 16 The fully localized limit (FLL) was adopted, 
as formulated by Liechtenstein et alJ^- In addition, we 
performed GGA+U calculations as well. In both cases, 
we neglect multipole components in the expansion of 
the Coulomb interaction by using only the U parameter, 
and setting the intra-atomic exchange parameter J = 0. 
Thus, U e s = U — J = U , and the Slater integrals F 2 and 
F A are zero. In LDA+C//GGA+C/ calculations, one of- 
ten chooses the U parameter such, that the experimental 
magnetic moment or band gap is reproduced. 

As another route to obtain an electronic structure with 
the full magnetic moment, we use the PBE functional 
and enforce the experimental (or Slater-Pauling) magne- 
tization within a fixed spin moment (FSM) calculation. 
Thus, we obtain the (variationally unique) PBE band 
structure with the lowest energy, that results in the de- 
sired magnetization^ This method is usually used to 
map the dependence of the total energy on the magneti- 
zation. Here, we use it to model the unoccupied states of 
Co2FeSi, such that the minority d-states above the gap 
are shifted above the Fermi energy (see the upper part 
of Figure [6] for a comparison of GGA and GGA FSM re- 
sults). We have performed the calculations with a fixed 
total spin moment of 6 /iB / f-u. 

The muffin-tin radii Rmt were set to 2.2 a. u. for all 
atoms, and the plane- waves cutoff was set to i?MT'^max = 
8.0. The angular momentum expansion in the muffin- 
tins was taken to Z max = 10 for both the wave functions 
and the potential. The maximum G- vector length for the 
density and potential expansion in the interstitial was set 
to |G| max = 14.0 a.u. -1 . k-point meshes of 20 x 20 x 20 
points, restricted to the irreducible wedge of the Bril- 
louin zone, were used. Spin-orbit coupling was included 
by the second-variational method. The 2p core states 
were described as (scalar relativistic) local orbitals, i.e., 
the 2p 3 /2 and 2pi/ 2 state have the same basis function. 
This can give rise to a small error for the 2pi/ 2 state, 
when compared to a fully relativistic description of the 
core states. Exchange splitting of these states is taken 
into account. X-ray absorption spectra were computed 
within the Kubo linear response formalism as lined out 
in Ref. [l9j. All theoretical spectra are broadened with a 
0.3 eV wide (FWHM) Lorentzian to account for lifetime 
broadening. 

Experimental lattice parameters of CFS were reported 
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TABLE II. Structural properties of the samples: out-of-plane 
lattice constants, B2 and L2i order parameters, and rocking 
curve full width at half-maximum. 



Name 


a (A) 


<Sb2 


Sl J 2 1 


Auj 


CFS-Cr 


5.62 


0.97 ±0.02 


0.8±0.05 


1.6 


CFS-HTD 


5.65 


0.98 ±0.02 


1.0±0.05 


0.5 



TABLE III. Element resolved magnetic spin and orbital mo- 
ments of CFS-Cr obtained from sum rule analysis of the TEY 
and luminescence spectra. All moments are given in /ib- 









mi 


mi/m s 


Co 


1.32 ±0.2 


1.38 ±0.2 


0.10 ±0.01 


0.077 


Fe 


3.17 ±0.4 


3.01 ±0.4 


0.09 ±0.01 


0.030 



to be 5.64 A and 5.66 A£ii We used the average value 
(5.65 A) for our calculations. 



IV. RESULTS 
A. Structural characterization 



All characteristic features that have been observed by 
other authors for Co2FeSi are present in our XAS/XMCD 
spectra as well: the shoulders on the trailing edges of the 
absorption edges, and the feature in the Co spectra 4eV 
above the absorption onset The latter arises from an 
s-d hybrid state of Co and Si. The luminescence spectra 
(not shown) are in good agreement with the TEY data. 

The element specific magnetic moments are extracted 
from a sum rule analysis of the spectra, assuming the 
radial p — >• d transition matrix constants C Fc = 6.6 eV 
and C Co = 7.8 eV according to Ref. [M The spin 
moment ratio determined from TEY is m^ c /m^° = 
2.4 ± 0.3, whereas the ratio from luminescence detection 
is m^/rrig = 2.2 ±0.3. Both samples have extrapolated 
low-temperature magnetizations of 6.1 ±0.2 /ie / f.u. (de- 
termined by VSM); using this value to scale the XMCD 
sum rule results, we find element-resolved moments as 
given in Table IIIII The scaling removes experimental un- 
certainties with respect to the magnetic saturation of the 
samples. 

With the white line sum rule from Ref. [22|, we obtain 
the numbers of d holes: 2Vj°° = 3.0 ± 0.4 and N% e = 
4.3 ± 0.6, which are significantly higher than the values 
for the pure elements (2.4 and 3.45, respectively)^ The 
analysis of the sum rules rests on the assumption that the 



The samples were epitaxial with a [001] orientation. 
The out-of-plane lattice parameters are close to the bulk 
values (Table |TT| . The structural ordering is given by the 
long-range order parameters Sb2 and Sl2i ■— 

L2i order is observed for both samples, and the order 
parameters Sb2 are close to 1 for CFS-Cr and CFS-HTD 
(Table HI)). The x-ray scattering amplitudes of Co and Fe 
are very similar, thus we can not detect chemical disorder 
between Co and Fe. Therefore, <Sb2 measures the disor- 
der between Si and Co sites, which is virtually zero. Sl,2i 
measures the disorder between Si and Fe sites; we observe 
a high degree of L2i order for both samples. While CFS- 
HTD has approximately full L2i order, we find a slightly 
reduced ordering for CFS-Cr, with about 10 % of the Fe 
sites occupied by Si. Note the larger experimental uncer- 
tainty for Sl2i due to defocussing of the x-ray beam in 
off-specular measurements. 

CFS-HTD has an x-ray reflectometry profile that in- 
dicates island growth. The real structure of this film is 
similar to what was observed by Schneider et al.^ i.e., 
a discontinous film. CFS-Cr has a smooth surface. The 
rocking curve widths Auj (reflecting the crystalline orien- 
tation distribution) are very different for the two samples, 
being small for CFS-HTD and rather large for CFS-Cr. 



B. X-ray magnetic circular dichroism 

In Fig. Q]we exemplarily show the x-ray magnetic cir- 
cular dichroism spectra obtained in total electron yield 
on sample CFS-Cr. 
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FIG. 1. X-ray absorption and magnetic circular dichroism 
spectra of the sample CFS-Cr taken in total electron yield. 
The absorption spectra are normalized to unity 40eV above 
the L3 absorption onset. 
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FIG. 2. Fundamental x-ray magnetic linear dichroism spectra 
of Co in Co2FeSi, taken in total electron yield. The normal- 
ization is the same as in Fig. [1] 
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FIG. 3. Fundamental x-ray magnetic linear dichroism spectra 
of Fe in Co2FeSi, taken in total electron yield. The normal- 
ization is the same as in Fig. [1] 



radial p — > d transition matrix constants are independent 
of the chemical environment. This might be true within 
about 10 % and gives rise to a corresponding error. Fur- 
ther, there is some uncertainty due to the background 
subtraction, which we estimate to be of the order of 10 % 
as well. 



C. X-ray magnetic linear dichroism 

The XMLD is generally measured as the difference be- 
tween two absorption spectra taken with linearly polar- 
ized light; one with the polarization vector parallel to the 



magnetization, and one perpendicular. Expressed with 

TEY intensities, i.e. XMLD = /tey~ 7 tey- The XMLD 
has a large crystal anisotropy. It was shown for cubic sys- 
tems that all XMLD spectra can be described as a linear 
combination of two fundamental spectra, which are mea- 
sured along the [100] and [110] high symmetry directions, 
respectively^ Earlier reports of XMLD of Heusler com- 
pounds include Co 2 MnSi and Co 2 MnAl, 25 Mn 2 CoGa and 
Mn 2 VGa^ and Co 2 TiSn^l 

The XMLD spectra of CFS-Cr are presented in Fig- 
ures [5] and The spectra have pronounced substruc- 
tures, and the intensities are quite large. For Co, the 
maximum XMLD intensity is found in the [110] direc- 
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tion at the L3 edge, and amounts to 7.1% with respect 
to the L 3 resonance height. The maximum intensity for 
Fe is observed at the L 3 edge in the [100] direction, and it 
amounts to 6.9%. For Co, the maximum XMLD is twice 
as large in the [110] direction than in the [100] direction. 
In contrast, for Fe the signal strengths are similar in the 
two fundamental directions. 

Comparing the Co XMLD spectra with the available 
published data of other Heusler compounds, we note that 
the [110] spectra are very similar in all cases (Co2MnSi, 
Co 2 MnAl, Co 2 TiSn, and Mn 2 CoGa)A-^ The only ex- 
ception is Co2MnSi, which has a sharper, and more de- 
tailed structure. Telling et al. ascribe this (in particular 
in comparison with Co2MnAl) to the half-metallic char- 
acter of Co2MnSi and the associated localized character 
of the d electrons^ This could hint at the lack of half- 
metallic character in Co2FeSi. 

It has been demonstrated that the XMLD of Heusler 
compounds depends on the square of the spin moment, 
modified by the degree of localization of the magnetic 
moment 2§r2& The maximum Co XMLD intensity along 
[110] at the L3 edge as a function of the square of the 
Co spin moment (4.1%//ig, using the TEY values for 
consistency) is weaker than for Co2MnAl (4.79%//i|), 
Co 2 MnSi (7.6%//4), or Co 2 TiSn (5.93%//z|). The 
heights of the L3 resonances with respect to the back- 
ground are comparable among these compounds, so this 
normalization is reasonable. This hints at not fully 
localized Co moments, and indicates the lack of half- 
metallicity as well. Still, Co2FeSi has the largest XMLD 
observed for Co in a Heusler compound due to its large 
magnetic moment. 

The Co [100] XMLD spectrum of Co 2 FeSi is different 
from the spectra of the other compounds mainly by a pos- 
itive contribution at the maximum of the L3 resonance, 
which is very weak in the other cases. 

XMLD spectra of Fe containing Heusler compounds 
are not available in the literature. However, our spec- 
tra are clearly different from XMLD spectra of pure 
Fe films grown on GaAs(001)i22 Notably, the Fe spec- 
tra of Co2FeSi are similar to the Mn XMLD spectra of 
Co 2 MnSi. 2 - The [110] spectra have a double-peak con- 
tribution at the low-energy side of the L3 edge and a 
negative contribution at the high-energy side. The L2 
edge is even more similar. For the [100] direction we ob- 
serve a similar resemblance, but Mn has a small positive 
contribution at the absorption onset, which is missing 
for Fe. The maximum XMLD intensity of Fe as defined 
above (0.69 %//4) is close to that of Mn in Co 2 MnSi 
(0.77 or Co 2 MnAl (0.68 and larger than 

in pure Fe (0.54%//x|) . 28 i 29 Fe and Mn occupy the same 
site of the Heusler structure and experience similar lo- 
cal crystal fields. The larger XMLD intensity of Fe in 
Co 2 FeSi compared to pure Fe indicates a more localized 
character of the magnetic moment. These findings em- 
phasize the high site specificity of the XMLD and its 
sensitivity to the local crystal field. 

In Figure|l]we compare the XMLD spectra of our sam- 
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FIG. 4. X-ray magnetic linear dichroism spectra of the sam- 
ples CFS-HTD (red traces) and CFS-Cr (blue traces). The 
spectra were taken in luminescence detection. Subtle differ- 
ences in the intensities of the various features are visible (see 
arrows). 

pies taken in luminescence detection mode. The intensi- 
ties are given as percentages of the L3 absorption peak. 
The main purpose of this comparison is to verify that 
the XMLD spectra are essentially independent of the de- 
position conditions and specific details of the films. The 
various spectra are very similar throughout the samples, 
and only small differences are seen for the intensities of 
individual features. These differences may be associated 
with the slightly lower degree of L2i ordering in the CFS- 
Cr sample. 

D. First principles calculations 

Four computational approaches are compared with ex- 
periment: plain GGA, GGA fixed spin moment (FSM), 
LDA+J7 (U = 4.5 eV), and GGA+J7 (U = 1.9 eV). For 
LDA+[/ we adopt the U value of References and H3. 
In our GGA+C7 approach, we can use a slightly lower 
U than in LDA+[7 to obtain the full magnetic moment. 
In the latter case, at least 2.5 eV are required; here, we 
adopt the smallest value that results in the full magnetic 
moment. This value was also used in Ref. within a 
GGA+U approach. The total magnetic moment neglect- 
ing spin-orbit coupling with the plain GGA (with the 
PBE functional) is 5.53 /iB, which is in good agreement 
with previously reported values.^ Figure [5] shows the the- 
oretical and experimental spectra for Co and Fe along the 
[100] and [110] directions. 

The general shape of the GGA spectra is quite similar 
to the experimental ones, but the width of the theoretical 
curves is considerably too small. The GGA FSM descrip- 
tion improves considerably the width and general shape 
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FIG. 5. X-ray magnetic linear dichroism spectra calculated with different approaches, as discussed in the text. The theoretical 
spectra are vertically offset against the experimental spectra. They are aligned in energy and scaled to the experimental data. 



of the spectra, and gives a very good overall resemblance 
of the experimental data of Co [100] and [110], and Fe 
[100]. 

In contrast, both the LDA+U and GGA+U spectra do 
not describe the experimental data. For Fe [100], the L 3 
spectra are inverted, and the Co spectra have spurious 
substructures. Further, the Co spectra are broader than 
the experimental data. The computed d band width is 
consequently too large. 

The Fe [110] spectra are very similar in the GGA, GGA 
FSM and LDA/GGA+C/ approaches. All of them have 
a negative contribution at the absorption onset, which 
is positive in the experiment. Some features at the high 
energy side of the L3 edges in the experimental spectra 
are not captured by any theoretical approach - these are 
marked with arrows in Fig. [5J In the Co spectra, these 
features are associated with the Co-Si s-d hybrid state. 
Since our one-electron theory fails to describe these fea- 
tures, we conclude they are associated with many-body 
effects. 

Electron-hole correlations, in particular the formation 
of bound excitons (which lead to the so-called 'multi- 
plet' splittings in the absorption spectra), are often ob- 
served in the L32 x-ray absorption spectra of 3d tran- 
sition metal atoms in insulators. However, these effects 



are virtually absent in the spectra of metals, and reduce 
to a change in the branching ratio away from the statis- 
tical 2:1 in the early 3d transition metals. Laskowski 
and Blaha have shown by solving the Bethe-Salpeter- 
Equation in an all-electron framework that the dielec- 
tric response of the system plays a crucial role for these 
effects^ The formation of the splittings is governed by 
the screened Coulomb interaction of the excited elec- 
tron with the core-hole. In a (half-)metal, this inter- 
action is exponentially screened out, so that only small 
changes of the spectra remain. In a recent paper we have 
demonstrated that this effect is of the order 0.3 eV for 
Co in Co2TiSn by means of core-hole calculations in the 
independent-particle-approximation^ whereas the un- 
occupied d band width is of the order 2 eV here. Thus, 
due to the metallic screening of the Coulomb interac- 
tion it is reasonable to assume that independent-electron 
calculations can give a good description of the x-ray ab- 
sorption spectra of Co2FeSi. 

From the comparison of the calculations, it becomes 
clear that GGA FSM provides the best overall matching 
of the experimental data, and can serve to model the un- 
occupied band structure of Co2FeSi. Based on the GGA 
FSM electronic structure we have computed the number 
of d holes of Fe and Co. This has to be done with care in 
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the FLAPW method, because these numbers depend on 
the choice of the muffin tin radii. Therefore we computed 
the hole numbers with different choices of muffin tin radii 
and selected the ones for which the partial magnetic mo- 
ments saturate - the relevant part of the d orbitals is then 
presumably entirely contained in the muffin-tins. These 
radii (2.2 a. u.) were used for all calculations through- 
out this work (muffin tin radii of 2.28 a. u. correspond 
to maximum space filling). The corresponding numbers 
of d holes are N%° = 2.9 and N% e = 4.1. Both are a 
bit smaller than the values extracted with the white line 
sum rule (see Section I IV Bp . However, they are larger 
than the theoretical hole numbers of the pure elements, 
which amount to 2.4 and 3.45, respectively^ Our val- 
ues are close to the numbers of unoccupied d states of 
the free Co and Fe atoms. This atomic-like configuration 
could be the reason for the failure of the PBE functional 
to obtain the correct magnetic moment. 

The calculated magnetic moments within the GGA 
FSM approach are = 1.49 hb, mf° = 0.043 and 
mf e = 2.98 /xb, mf c — 0.043 /zb- The interstitial contri- 
bution is almost zero. The calculated spin moment ratio 
is 2, which is smaller than the experimental value. The 
orbital moments are underestimated by a factor of about 
2.3, which is typically observed in DFT calculations with- 
out orbital corrections^ 

We note that it is still an open question why bulk 
magnetometry of presumably half-metallic Heusler com- 
pounds gives the integer Slater-Pauling values, which are 
only valid in the complete abscence of valence-band spin- 
orbit coupling. XMCD does, however, always give small 
orbital moments, which should increase the total mo- 
ment. 



V. DISCUSSION 

Kunes et al. proposed a theoretical model of the 
XMLD, in which the XMLD is given as an energy deriva- 
tive of differences of the unoccupied e g and t2 g densities 
of states: 

d 



XMLD oc A x — [a{t2ff\ 
ail/ 



hgi) + 0(egx - e gi )] , (1) 



where A x is the core-level exchange splitting, and a, 
are constants, which take the values a = — 1 and = 2 
for the [100] magnetization direction! 29 ! 32 The energy de- 
pendence has been omitted for clarity. This model is only 
valid without spin-orbit coupling, which lifts the degen- 
eracies of the crystal-field split doubly and triply degen- 
erate d states. However, this model gives an instructive 
picture of the local crystal field sensitivity of the XMLD. 

Consequently, the large differences observed for the 
GGA FSM and GGA+C/ spectra along the [100] direc- 
tion indicate significant changes in the distribution of the 
orbital symmetries (irreducible representations) of the d 
states. 

In Figure [5] we compare the total DOS of the GGA 
FSM calculation and the GGA+U calculation neglecting 
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FIG. 6. Top: Total densities of states of Co2FeSi with the 
GGA fixed spin moment (FSM) and GGA+<7 (1.9 eV) ap- 
proaches. The insets show the vicinity of the Fermi level. 
The light blue line in the top panel represents the GGA cal- 
culation. Bottom: Partial densities of states of Co and Fe in 
the two computational approaches. The spectra are decom- 
posed into the irreducible representations (e g and t2 g ) of the 
d states. Spin-orbit coupling is neglected. 



spin-orbit coupling. The gap widths and positions of the 
Fermi levels are different, and the width of the unoccu- 
pied minority d states is larger in the GGA+J7 calcula- 
tion. The size of the minority gap in the GGA+f/ cal- 
culation is about {7/2. The GGA FSM density of states 
is essentially the same as that of a plain GGA calcula- 
tion, with the majority and minority states being rigidly 
shifted against each other. Due to their overall similar- 
ity, typical experimental techniques that probe the DOS 
(photoemission, inverse photoemission, x-ray absorption) 
are hardly able to distinguish between the models of the 
electronic structure. 

By decomposing the individual d projections into their 
orbital symmetries, we are able to observe differences 
between the GGA FSM and GGA+<7 models. This is 
shown in the bottom part of Figure [B] The DOS of the 
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two orbital symmetries (e g and t2 g ) have entirely differ- 
ent structures in the two calculations. The e g and t2 g 
states are more split in GGA+U (i.e., the crystal field is 
stronger), and the hybridization between Co and Fe t2 g 
states becomes stronger as well. The Co e g states are 
broader in GGA+U, and the Fe t2 g states have a larger 
weight in GGA+U than in GGA FSM. The crystal field 
splitting is the driving mechanism for the gap creation 
in the half-metallic Heusler compounds. 1 Obviously, the 
+U operator (in the FLL formulation) enhances the crys- 
tal field splitting and forms the half-metallic gap. 

It has been shown that the LDA+J7 formalism can 
be derived as an approximation to the quasiparticle self- 
energy within the GW method i 33 ' 34 The derivation in- 
volves three assumptions: (1) the screened Coulomb in- 
teraction (described by U) is static; (2) only localized 
states are subject to the quasiparticle corrections; (3) 
the many-body exchange interaction between localized 
and itinerant electrons is neglected. Thus, LDA+[7 is 
a kind of on-site Hartee-Fock theory with the effective 
screening given by the +U parameter. It does pro- 
vide results that are quite similar to those from the hy- 
brid PBEO functional^ The inclusion of Hartree-Fock 
exchange without the proper spacial and frequency de- 
pendence of the screening (i.e., the correlation) allows to 
obtain the magnetic moment of Co2FeSi correctly, but 
fails at the correct prediction of the unoccupied states 
band structure. 

The GGA is able to describe the wave functions and 
unoccupied states reasonably well up to a rigid shift of 
the latter. Since the GGA can be seen as a static, local, 
and orbital-independent approximation to the GW self- 
energy, the success of reproducing the XMLD spectra 
well is based on a cancellation of errors, due to omis- 
sion of both the non-locality and the frequency depen- 
dence in the exchange and correlation contributions. The 
good balance between exchange and correlation in the 
LDA (and GGA, which is essentially the same theory) 
has been confirmed earlier by Kotani and Griming et 
al. These authors have performed calculations of transi- 
tion metal and semiconductor band structures with the 
EXX+RPA approach, which treats exchange and cor- 
relation on the same level of many-body perturbation 
theory^^ EXX+RPA reproduced the LDA results well, 
including magnetic moments and exchange splittings of 
Fe, Co, and Ni. However, for MnO, where screening 
is weak, EXX+RPA gave substantial improvement of 
the band structured Furthermore, GW calculations on 
half-metallic CrAs and MnAs by Damewood and Fong 37 
demonstrated that the GGA can predict half-metallic 
gaps that are very close to the GW values, depending on 
the details of the local fields and the associated screening. 



The GGA FSM band structure has a small gap (about 
0.1 eV wide) in the minority states, and the Fermi energy 
is located in the gap. Inclusion of spin-orbit coupling 
closes this small gap. This observation may explain the 
rather low TMR values of magnetic tunnel junctions with 
Co2FeSi electrodes, the lack of half-metallic characteris- 
tics in the tunneling conductance curves, and the large 
Gilbert damping^ 

VI. CONCLUSIONS 

Our experimental results show that the features associ- 
ated with the GGA+U band structure are not present in 
Co2FeSi. GGA+U predicts an inadequate distribution of 
the crystal-field split d states, which we identified by x- 
ray magnetic linear dichroism. The +U operator opens a 
minority gap at the Fermi energy, which necessarily leads 
to an integer magnetic moment. While exchange and cor- 
relation effects beyond the GGA are certainly important 
for the description of Co2FeSi, the simple +U approach 
does not produce the correct band structure. The PBE 
functional is able to produce an unoccupied states band 
structure for Co2FeSi which is approximately correct up 
to a rigid shift of the minority d states. 

Co2FeSi is a complex material, being neither close to 
the electron gas (the GGA fails to find the correct ground 
state magnetic moment), nor being strongly correlated 
(GGA+U fails at the description of spectral properties). 
That calls for further refinement of the approximations 
to the density functional. More flexible functional forms, 
such as meta-GGA functionals (adding information on 
the kinetic energy density), may perform better in this 
respect. Co2FeSi is therefore an interesting test case 
for new density functionals. Ultimately, many-body ap- 
proaches, such as EXX+RPA, 35 the GW approxima- 
tion (possibly including some kind of self-consistency), 
or DMF T 38 i 39 might be required to obtain a consistent 
theoretical picture of Co2FeSi. 
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